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Abstract. The method of calculation of the resonance characteristics is developed for 



the metastable states of the Coulomb three-body (CTB) system with two disintegration 
channels. It is based on the numerical solution of the scattering problem in the framework 
of the adiabatic hyperspherical (AHS) approach. The energy dependence of .fí-matrix in 
the resonance region is calculated with the use of the stabilization method. Resonance 
position Eq, partial widths Ti, and three additional parameters are obtained by fitting 
of the numerically calculated Kij(E) with the help of the generalized Breit-Wigner for- 
mula which takes into account the non-zero background inelastic scattering. The method 
\ developed is applied to the calculation of the parameters of three lowest metastable states 

of the mesic molecular ion dtu. 

m 
o 

P Introduction 



The goal of the present paper is the extension onto the two-channel case of the numerical 



method developed previously for the resonances in the CTB systems with one open channel 
PP, 0. In particular, the special attention is devoted to the proper calculation of the 
q-( partial widths Tj which are interesting for applications together with the total width 

r = £rv 

The method developed is applied to three lowest resonances in the dtfi system with 
zero angular momentum (J = 0) and two disintegration channels: (1) (í/x) n =i + d and (2) 
(d/i) n= i + t. They represent the quasistationary states in the potential well formed by the 
third AHS term (fig.l). 

The resonances pointed out are interesting for many mesic atòmic processes and were 
investigated repeatedly. The total widths r = F 1 + T 2 were calculated in papers [3] jl] 
(complex rotation method), [ü] (analysis of elàstic and inelastic scattering), [ü] (Siegert 
pseudo-states method), [7j (i?-matrix method). The results of [3] jl] differ significantly 
(3-r5 orders) from those of jH], jüj, [Zj- The results for the partial widths which play an 
important role in the theory of the muon-catalyzed fusion were presented up to now only 
in paper [5] where the non-adiabatic coupled rearrangement channel method S] had been 
used for the calculation of the elàstic and inelastic cross sections. 

Our method is based on the AHS approach which is widely used in the calculations 
of mesic atòmic systems starting from paper |ÏÏ] (see for references (Ull)- The three-body 
wave function is represented as a series over AHS basis, and the determination of the 
reaction matrix is reduced to the numerical solution of the system of radial equations 
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Figure 1: AHS terms £j(p) for the dtfi molecule. The resonances under consideration are 
connected with the potential well in the third term. 

(sec.l). It is calculated as a function of the energy E in the resonance region by the 
methods described in paper pQ. In particular, the stabilization method [11] is used. 

The resonance parameters are obtained by the fitting numerical results for Kij(E) 
with the help of the extended Breit-Wigner formula presented in section 2. It takes into 
account the nonzero background inelastic scattering and therefore contains (n+l)(n+2)/2 
independent parameters (n is the number of open channels) instead of (2n + l) parameters 
in the well-known traditional expression Moreover, this formula imposes the 

non-trivial restrictions on the vàlues of IV As it is seen from the results for the lowest 
resonance in dt/i presented sec. 4 the use of the extended formula in this case is necessary 
for an adequate representation of Kij(E) and for a proper calculation of partial widths. 

The method developed relates to ones based on the numerical solution of the scatter- 
ing problem. The typical feature of these methods is the possibility to obtain the high 
accuracy for T while E is calculated not so precisely p. Indeed, in the scattering prob- 
lem E and T appear as the vàlues of the entirely different nature: E is the position of 
the resonance peak while T _1 is its height, the errors of T and E are the independent 
vàlues in a sense. On the other hand in the complex rotation method the resonance po- 
sition Eq and the half-width T/2 appear as the Cartesian coordinates of the pole in the 
complex .E-plane, they are considered usually as the vàlues having the same rights and 
are calculated with the same absolute error. 
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1 Determination of ií-matrix in AHSA 

We use the mesic atòmic units (h = e = m M = 1) and Jacobi coordinates 

m t r t + m d r d 

R = r d -v t , r = ■ . (1) 

m t + m d 

The hyperspherical coordinates p G [0, oo), x £ [O 71 "]) $ e [0, tt] are defined by formulae 

p 2 = 2pr 2 + 2MR 2 , 

, X //í\ 1/2 r q rR 

tan — = — — , cosit = — — , 
2 \MJ R 1 rR 1 

fT 1 = 1 + (mt + ma)- 1 , M~ l = mj 1 + . (2) 
The three-body Hamiltonian H in the case J = has the form jÜ] 

where the adiabatic Hamiltonian /z is given by the expression 

4/<9 d 1 d d \ 

h = ^-^( — sm 2 X ^ + —^ 1 ^smiï—)+V(p,x,iï), (4) 

p 2 sm x V c/x ^X sm v c/ï/ c/ï/ / 

and V is the sum of the Coulomb interactions between t, d and /1. 

The three-body wave function ^(R, r) in the case J = depends only on p, x, 
It is presented in the form of the AHS expansion 

00 

r) = p~ 5/2 ^2 fj(p) ( Pj(p\Xj (5) 
i=i 

where the AHS basis functions ^-(p|x, 1?), j = 1,2, ...,00 are the eigenfunctions of the 
AHS Hamiltonian h (4): 

hfj{p\x, = £j{p)<Pj{p\x, (6) 
Their normalization is defined by the relation 

7T 7T 

(^iki) = / J í Pi(p\x,^) í Pj(p\x,^) sm 2 xsmddxd'd = 5ij. (7) 

o o 

Three lowest AHS energy terms £j(p), (j = 1, 2, 3, 4) are presented in fig.l, where their 
limiting vàlues at p — > 00, i.e. the energies of bound states of corresponding àtoms, are 
pointed out: 

ei )3 (oo) = E^in = 1, 2), £ 2i4 (oo) = E dlx {n = 1, 2). (8) 

Radial functions fj(E\p) satisfy the infmite system of coupled ordinary differential equa- 
tions 
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+ E 



d d 

Hjjifj, + Qjy — Jy + — (Qjj'fj' 



(9) 



(10) 



To calculate the two-channel reaction matrix Kij(E) = Kji(E), = 1,2) corre- 
sponding to channels (1) tp + d and (2) dp + t with thresholds E = £i(oo) and E = 62(00) 
one has to obtain two linearly independent solutions of system (9) (we label them with 
upper index i = 1,2) with boundary conditions (qf = E — e^oo)) [Tüj : 



fí(p) 



P 



0, j = 1,2, ...,00; 



/;(p) 



0, j = n + 1, n + 2, 00; 



" J " ' p — > 00 

% sin(çjp - ttJ/2) + (qi/qjf^Kij cos(q ó p -nJ/2), j 



f}(p) ^ ^Sa àiHu.P - 7T-//2) + (</,•/</, I 1 - "A';,- cosí./,/; - ttJ/2). / 1.2. (11; 

The partial cross sections of elàstic (i = j) and inelastic (i 7^ j) s-scattering (J = 0) 
is expressed in terms of (left idex corresponds to input channel) jTÜ] : 



4tt 5 íj D 2 + K 2 



k? (1 - D) 



F 2 



D = K 11 Kn-K 12 K 2h 



i,j = l,2, 



F = K U + K 22 , 



(12) 



fei = (2/i,) 1 / 2 ^ 



(mt + 1) 1 + 



ír'+m: 1 . 



2 Generalized Breit-Wigner formula 

The general formula describing the n-channel scattering matrix as a function of energy E 
in the vicinity of isolated complex pole Eq — iT/2 with small imaginary part (T/Eq <C 1) 
can be derived from the unitary property, the symmetry and the supposition that the 
residue of S'-matrix in this pole has an order 0(T). In the special case when the non- 
resonant inelastic scattering is absent and the background term has a diagonal form (we 
use the symbols with tilde for this special case) this formula (the Breit-Wigner formula) 
reads [□], 03] : 

S(E) = S b -- ' — , (13) 

v ; e-e q + íy/2' v ; 

where the background S'-matrix S b and the residue matrix B do not depend on E and 
are defined by the expressions 

S>. = 8 ió e 2iA % B l3 = (14) 

Here f3{ (i = 1,2, ...,n) are the real parameters (positive and negative) saisfying the 
condition 

n 

EÀ 2 = i- (15) 

1=1 
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Thus, in the absence of the background inelastic scattering the Breit-Wigner formula 
contains (2n + 1) independent real parameters: E , T, independent phaseshifts Aj (j = 
1,2, ...,n), and parameters $j (j = 1,2, ...,n) connected by relation (13). The value 

n 

r^rpl j = i,...,n, (Er, = r) (ie) 

is the partial width corresponding to channel j. 

In the general case, when the background inelastic amplitude is not a negligibly small 
value, it is necessary to modify formula (13). One has to find the general expression for 
matrix B in the formula 

S(E)=S b — -, (17) 

where S b is the symmetric unitary matrix which may be the nondiagonal one. This 
derivation can be easily done with the use of the orthogonal matrix R reducing S b to the 
diagonal form: 

S b = R-^t'R. (18) 

Here the diagonal matrix S b is given by the first of eqs.(12), so {k = 1, ...,n) are the 
eigenphases of matrix S b . Now the orthogonal transformation S(E) = R~ 1 S(E)R reduces 
the problem of the search for the general form of S(E) to the special case (13), and as a 
result we obtain the general expression for S(E) in the resonance range in the form: 



S(E) = S b — = R" 1 

K ' E-E + iT/2 



E - E + iT/2 



R, (19) 



where S b and B are given by eqs. (14), (15). 

The general expression for S(E) (19) thus contains (n + l)(n + 2)/2 real independent 
parameters: E , T, independent phaseshifts Aj, parameters f3j connected by eq.(15), and 
n(n — l)/2 parameters of an orthogonal matrix R. 

The matrix elements of B in eqs. (17), (19) have the form 

n n 

B jk = E RijBi m R mk = E RijPie lAl E R mk P m e tAm = hfae^ +Ak \ (20) 

l,m 1=1 m=l 

where f3j (j = 1, n) are the complex vàlues connected with parameters (3 by the relation 

n n 

^• = EW (A - A ^, (EIAl 2 = i). (21) 
i=i i=i 

The probability of the disintegration to the channel i (partial width Tj) is equal to 

n 

r í = r|AI 2 , i = i,...,n, (E r ^ r )- (22) 

í=i 

It coincides with "eigenwidth" Tj only if R is an identity matrix. 
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The essential difference between the set {Tj} and the set {Tj} is that at fixed back- 
ground scattering matrix and complex pole (i.e. at given R, Aj, E and T) one can to 
define the positive parameters r\ arbitrary (with only restriction Y%=i = T), while the 
domain of definition of the parameters Tj depends on R. For example, in the case n = 2 
the transformation matrix has the form 



R = 



( cos v sin v \ 



\ — sin v cos v ) 



(23) 



where v is the mixing parameter. The analysis of the expression for $ (21) in this case 
leads to the following restriction for Ti/T = \[3i\ 2 : 



2iyr > 1 - a/cos 2 (Ai - A 2 ) sin 2 2v + cos 2 2u > 0. 



(24) 



This inequality means that in the presence of the background inelastic scattering [y ^ 0) 
the probability of the disintegration certainly differs from zero for every open channel. 

The complex pole (E — iT/2) of S- matrix corresponds to the real pole E 1 of the 
reaction matrix 

K = i(l + S)- 1 (l-S) = K b - ¥ Z^- ( 25 ) 

It is easy to obtain the expressions of E\, background reaction matrix K b , and matrix C 
in terms of parameters E , T, (3j, Aj and transformation matrix R: 



K b = R~ 1 K b R, 
C = R ^CR, C 



j'=i 

r PA 

2 cos Aj cos A k 



(26) 

(27) 
(28) 



If we have to calculate E , Tj, T using the numerically obtained K(E) we have to fit 
K(E) with the help of (n+ l)(n + 2)/2 real parameters in accordance with eqs.(25)-(28). 
In the case n = 2 it is necessary to determine six such parameters. We can use also 
six alternative parameters, namely, the pole of ií-matrix E\ and the independent matrix 
elements K\ x = a±, K\ 2 = a, K\ 2 = a 2 , Cu = bi, C 2 2 = b 2 : 



K 




(E - E,) 



ò \ 

&2 / 



(29) 



where 6i > 0, b 2 > 0, bib 2 - ò 2 = 0. 

The case of the zero background inelastic scattering corresponds to a — 0. 

Resonance position £o and widths Tj, T can be expressed in the explicit form in terms 
of E u a h bi with the help of eqs.(14)-(16), (26)-(28): 



Eq — Ei — 



h(l - d) 



d) 2 + f 



(30) 
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J=0 



A.(a) 




Figure 2: Eigenvalues of the auxiliary boundary problem Aj for J = as a function of 
right boundary a. The avoided crossings at Aj = E appear as exact crossings on the 
scale of the figure. 



li 



hü2 + bi — dbi 



ha\ + b 2 — dbi 



d) 2 + p ' 1 (1 - d) 2 + f 2 
Jh + g(l-d) 



Here 



r = r! + r 2 



d = a±a2 — a z 



;i - dy + p 



f = ai + a 2 , 



g = b\ + 6 2 , /i = aib 2 + a 2 bi — 2ab. 



(31) 
(32) 

(33) 



3 Numerical method 

Basis functions (fj(p\x, $) (6), (8) and matrix elements (10) were calculated in accordance 
with algorithms described in papers jH], ^3]- The calculations were performed on the 
orthogonal finite-element grid [A^ x x N#] with using the second order Lagrange elements. 
The number of nodes in x and $ was taken equal to N x = 131 and N$ = 61. This provided 
the accuracy of calculation ~ 10~ 5 of all matrix elements. 

The final results have been obtained with the number = 6 of AHS basis functions. 

The energy dependence of K- matrix (11) in the resonance range was obtained with 
the use of stabilization method jTTj in the same way as in papers 0, 0- We investigate 
the discrete spectrum of the auxiliary eigenvalue problem for radial system (9) on the 
finite interval < p < a; its eigenvalues Aj(a) have the avoided crossings at Aj(a) w i?o 
(see fig.2); we solve the scattering problem for radial system (9) and calculate K(E) at 
E = A.j(a) in the neighbourhood of the resonance; the scanning along a instead of along 
E allows to study in details the behavior of K(E) near E = E (see fig.3). 
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(E-E )*10 9 [m.a.u] 



Figure 3: Matrix elements of K-matúx as functions of E (J = 0, v = 0). 

The system of radial equations (9) was solved at p G [0.05, p* = 500 with adapted step 
|14j . Tiris provides tlre results with the relative accuracy ~ 10~ 4 . 

Parameters Ei, O-x, CÍ2, O, &lj 62 ; b of ií-matrix (11), (29) were obtained by fitting 
numerical results for K(E) with the help of the generalized Breit-Wigner formula (sec. 2). 
The relation b\b^ 2 = 1 was vàlid in this case with accuracy 10~ 4 . Parameters E , T, F\, 
and r 2 were calculated by formulae (30)- (33). 

As a whole we estimate the relative accuracy of calculation of T and Tj as 20 -4- 20%. 

4 Results and discussion 

The main results for three lowest resonances ( J = 0, v — 1, 2, 3) are presented in figs. 3-5 
and in tables 1, 2. 

The numerically calculated matrix element K\\ (v = 0) as function of (E — Eq) is 
presented in fig. 3. Dots correspond to the scanning with fixed step along a (sec. 3). The 
step along E in this case becomes very small near resonance. The use of the scanning 
along a allows to trace the .E-dependence of K\\ in details. 

Fig. 4 presents the numerical data for [K\\ — cti) -1 , {Kvi — a)" 1 , and (JÍ22 — «2) -1 as 
functions of {E — Eo) for v = 0. It is seen, in particular, that all matrix elements of 
K- matrix have the pole in the same point E . 

The numerically obtained elàstic cr n , cr 2 2 an d inelastic <j\i cross sections (12) are 
presented on fig. 5 as functions of the energy. All curves have a typical Fano form. The 
parameters of these three profiles can be expressed in terms of E\, a\, ci2, a, 61, and 62- 
The numerical data for <j\i show clearly that it can not be presented without taking into 
account the nonzero background scattering. Indeed, at a = the profile crui^E) should 
be the symmetric one with respect to E = E . 

The data presented in tables 1 show that for three lowest resonances in the dtp system 
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(E-E/10 9 [m.a.u] 



Figure 4: The inverse vàlues of matrix elements of the resonant part of ií-matrix ( J = 0, 
v = 0). 

the use of N = 6 AHS basis functions gives for E Q the accuracy 0.01%. To improve the 
accuracy and to perform the calculations for large J and v it is necessary to increase N 
and p*. 

Our results for T (table 2) coincide with those of papers [B], [Z| with an accuracy 
10-i-50%. Independently of the precision of the numerical calculation of reaction matrix 
K(E) the results for branching ratio can change essentially (about 10% in our case) if 
one uses for data processing the Breit-Wigner formula without taking into account an 
inelastic background scattering, i.e. if one assumes a = in eqs.(3)-(33). 
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Figure 5: Elàstic and inelastic cross sections in the resonance range (J = 0, v — 0). 

Table 1. Resonance position —E relative to n = 2 level of tfi (in /za.u.) 
for lowest resonances (J = 0,t> = 0,l,2)in dt/j, . 



V 


□ 


m 





Present 






N=40 


N=30l·l-50 


N=6 





0.159194 


0.1591938 


0.1591939 


0.15917 


1 


0.145302 


0.1453015 


0.1453019 


0.14524 


2 


0.134530 


0.1345291 


0.134529 


0.13445 



Table 2. Total width T (10 9 /za.u.) and the ratio r 2 /T (in parenthesis) 
for lowest resonances ( J = 0, v = 0,1,2) in dt\i 



V 





m 


[3 




□ 


Present 










N=40 


N=30^50 


N=6 





0.64xl0 4 


0.36xl0 3 


0.640 (0.1) 


0.354 


0.341 


0.47 (0.22) 


1 


0.14xl0 5 


0.50xl0 4 


5.56 (0.9) 


0.839 


0.829 


0.53 (0.38) 


2 


0.20xl0 6 


0.12xl0 5 


37 (0.9) 


1.15 


1.14 


1.02 (0.55) 
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